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The low temperature excitations in the anisotropic antiferromagnetic Fei_ x Zn x F2 for x = 0.25 
and 0.31, at and just above the magnetic percolation threshold concentration x p = 0.25, were 
measured using inelastic neutron scattering. The excitations were simulated for x — 0.31 using a 
localized, classical excitation model, which accounts well for the energies and relative intensities of 
the excitations observed in the scattering experiments. 

PACS numbers: 75.40.Mg, 75.50.Ee, 78.70.Nx 



I. INTRODUCTION 

In two and three dimensions, spin wave excitations 
are well studied in pure isotropic and anisotropic insu- 
lating antiferromagnetJ^l Magnetic excitations are sig- 
nificantly modified by magnetic dilution introduced by 
site substitution of the magnetic ions with diamagnetic 
ones. For isotropic systems near the magnetic perco- 
lation threshold concentration, x p , both well-resolved, 
local spin excitations as well as crossover from spin 
wave excitations to fracton excitations on a fractal-like 
lattic e 6 " 11 !, have been characterized near x p in two and 
three dimensions. Local spin excitations have been ob- 
served in dilute two-dimensional anisotropic systems near 
x p ^. In three dimensions, magnetic excitations have 
been studied for the magnetically dilute anisotropic sys- 
tems Mno. 5 Zn . 5 F^and Fe x Z ni _ x FjSHHl. 

The parent compounds MnF 2 and FeF 2 exhibit compa- 
rable exchange energies and corresponding spin wave dis- 
persions, but the FeF 2 system has an order of magnitude 
larger anisotropy and a correspondingly larger spin wave 
gap. The excitations in the Mno.5Zno.5F2 system were 
interpreted in terms of spin wave to fracton crossover as 
the scattering wavevector q increases. The behavior of 
the Fe x Zn!_ x F 2 , for x > 0.31, has been interpret ecP^ 
as showing both spin wave and local spin excitations for 
small q and local spin excitations for large q. In this 
study, we examine magnetic excitations with high reso- 
lution neutron scattering experiments and computer sim- 
ulations in Fe x Zni_ x F 2 as x approaches x p = 0.25. 

The well-characterized random-exchange antiferro- 
magnet Fe x Zni_ x F 2 , with its simple structure and inter- 
actions, is an ideal anisotropic three-dimensional (d = 3) 
system in which to study magnetic excitations through 



inelastic neutron scattering measurements and theoret- 
ical modeling and simulation. Magnetic excitations in 
the d = 3 anisotropic antiferromagnet FeF 2 have been 
very well characterized.^^ The structure of FeF 2 and 
diamagnetic ZnF 2 are similar.^ The antiferromagnetic 
spins in FeF 2 form a tetragonal lattice with two in- 
terpenetrating sublattices. The dominant antiferromag- 
netic inter-sublattice J 2 exchange interaction is between 
the body-center and body-corner magnetic ions. The 
intra-sublattice ferromagnetic J\ and frustrating antifer- 
romagnetic J3 exchange interactions are much smaller 
(cf. Sec. Ill for further details). Best fit values from 
inelastic neutron scattering measurements are shown in 
Table I. FeF 2 and diamagnetic ZnF 2 mix well during 
crystal growth to form Fe x Zn!_ x F 2 , which is a dilute, 
anisotropic, three-dimensional antiferromagnet. The oc- 
cupation of sites by Fe 2+ ions with S = 2 or diamagnetic 
Zn 2+ ions appears close to random, though slight clus- 
tering cannot be ruled out. It appea rs that J 2 does not 
vary significantly with dilution! 19 ! 20 ! There is limited in- 
formation about the effect of dilution on the anisotropy, 
but it also does not appear to vary by a large amount P3 

Magnetic ordering in F e x Zni _ x F 2 has been experimen- 
tally studied previously^ M 16 * 21 ^^! at magnetic concentra- 
tions x equal to or near x p — 0.246, the magnetic percola- 
tion threshold for the body-centered tetragonal magnetic 
structure with an interaction between the body-centered 
and corner ions (J 2 in Fe x Zni_ x F 2 ). The H = random- 
exchange transition should be expected for x > x p ii 
there is only the dominant exchange interaction J 2 . How- 
ever, the small J\ and J3 interactions in Fe x Zni_ x F 2 
could become influential near x p . The prior experiments 
in zero-field have demonstrated that for concentrations 
x < 0.31^ there is, at best, very weak long-range an- 
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tiferromagnetic order at low temperatures. The system 
exhibits spin-glass-like behavior, dominated by slow dy- 
namics near the percolation threshold, possibly a result 
of the frustrating J3 interaction. 

Early inelastic neutron scattering measurements in 
Fe x Zni_ x F2 were compared to a simple treatment with 
the excitation energies assigned^ as (z/n)E(q), where 
n = 8 represents the number of neighbors in FeF 2 , z < n 
is the possible number of neighbors of a given spin in 
the magnetically dilute system, and E(q) is the spin 
wave energy as a function of the scattering wave vec- 
tor q in Fe x Zn!_ x F 2 . The intensities are assigned by 
the combinatorial probabilities of finding z neighbors of 
a given spin. While giving a fairly accurate description 
of the overall spread in energy, this description fails in 
the detailed structure of the excitation spectrum when 
higher energy resolution measurements resolve individ- 
ual peaks. Similar results were found in far-infrared ab- 
sorption experiments, high magnetic field pulsed laser ab- 
sorption, and inelastic neutron scattering experiments for 
x > 0.4p3H31] it was observed in the pulsed laser absorp- 
tion measurements that the peaks become more easily 
resolvable for x = 0.4 and that for this case the simplis- 
tic modeling described above proves wholly inadequate^ 
the spacings of the resolved peaks do not correspond 
to the simple model. The excitations were found to be 
largely localized, having little dispersion. 

Here, we present a high resolution neutron-scattering 
study of Fe x Zni_ x F2 close to its percolation threshold. 
In agreement with the aforementioned work, the exci- 
tations show little or no dispersion. We show that a 
model of localized excitations accounts for our spectra 
even quantitatively. 

The layout of the remaining part of this paper is as fol- 
lows. In Sec. [TlJ we report our new experimental results. 
The obtained spectra are rationalized through a simple 
model in Sec. |III| Finally, we present our conclusions in 
Sec. El 



II. EXPERIMENTAL RESULTS 
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FIG. 1: Experimental intensities vs energy for x — 0.25 and 
0.31 for the zone center, q — 0, for q = 0.25, and for the 
zone boundary, q — 0.5. The data for x = 0.25 were taken 
at T = 0.7K and those at x = 0.31 were taken at T = 10K. 
The curves are guides to the eye constructed from Gaussian 
peaks. The Gaussian peak locations are the same for each 
value of q. The ratio of intensities for x — 0.25 and x — 0.31 
is arbitrary. 



tals were determined using density measurements. The 
resulting scattering spectra of the experiments, shown 
in Figs. [T] and [2j indicate localized excitations since the 
excitation energies are largely independent of the scatter- 
ing wave vector q. The simple approximation described 
above yields peaks similar to the resolved experimental 
peaks in the spectra, but the energies are not well pre- 
dicted, as discussed below. It was clear that a more re- 
alistic calculation was needed to describe the peaks and 
to elucidate what governs the details of the E vs q spec- 
tra. We discuss simulations below that capture essential 
characteristics of the experimental results. 



The results at x = 0.4 motivated experiments closer to 
the percolation threshold and we have conducted inelas- 
tic neutron scattering studies for x = 0.25 and x = 0.31. 
Neutron scattering measurements were carried out us- 
ing the high energy-resolution triple-axis spectrometer 
Cl-1 installed at the JRR-3M reactor of JAEA in Tokai 
operating with a horizontally focusing analyzer with a 
final neutron energy of Et = 3.1 meV. The energy res- 
olution at the elastic position is 0.09 meV (full width 
at half maximum) but it increases to 0.46 meV as the 
energy transfer increases to 8 meV. Single-crystal sam- 
ples were mounted in a closed-cycle refrigerator with the 
c-axis perpendicular to the scattering plane. We exam- 
ined a Feo.25Zno.75F2 single crystal with a mass of 2.24g 
and a Feo.31Zno.69F2 single crystal with a mass of 1.72g. 
The magnetic concentrations of the optical-quality crys- 



III. MODELING THE SPECTRA 

We have remarked in the previous section that, to a 
large extent, the scattering spectra are independent of 
the scattering wave vector, which suggests that the un- 
derlying excitations are spatially localized and therefore 
can be described by a local model. Local spin Ha miltoni- 
ans for pure FeF2 have been extensively discussed) 1 * 17 * 30 ! 
and include both an on-site interaction characterized by 
an anisotropy parameter D, and Heisenberg exchange in- 
teractions characterized by three coupling constants J\, 
J 2 and J3. The strongest exchange is the antiferromag- 
netic J2 that couples nearest neighbors belonging to dif- 
ferent sublattices of the body-centered tetragonal mag- 
netic lattice of Fe atoms, while Jj and J3 couple atoms 
belonging to the same sublattice, i.e., bonds parallel to 
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FIG. 2: Experimental intensities vs energy for x = 0.25 at 
T = 0.7, 20, and 40K to show the temperature dependence of 
the neutron scattering spectra. 





TABLE I: Best fit parameters, taken from Ref. 1, for the spin 
Hamiltonian of pure FeF2. 



Z^cm' 1 ) 



Ji(cm 1 ) J2(cm 



J3(cm 



6.46(+0.29,-0.10) -0.048 ±0.060 3.64 ±0.10 0.194 ±0.060 



sites in Fig. [3|. 

Tabic [I] shows the parameters obtained by Hutchings et 
aZ.E by fitting inelastic neutron scattering data of FeF 2 
to the spin Hamiltonian, and it is often assumed that the 
same values can be used for the analysis of the dilute an- 
tiferromagnelp3E^ (we will elaborate on this point later) . 
The fact that | J\ | , J3 <C J 2 < D allows us to get an 
estimate of the spectrum at low temperatures by ignor- 
ing in the Hamiltonian (TH) the terms proportional to J\ 

(z) (z) 

and J3 and making the approximation Si ■ Sj ~ Sj 
in the terms proportional to J 2 . For the sake of defi- 
niteness let us consider a site i on the A sublattice, with 
n 2 Fe neighbors (the probability distribution function for 
n 2 = 0, 1, . . . , 8 is binomial). Since at very low tempera- 
tures the magnetic state of the sample is essentially the 

( z) 

Neel state (S\ = ±2 if the site i belongs to the A sub- 
lattice and = — 2 if i is in the B sublattice), the local 
magnetic field felt by the spin i due to the n 2 surround- 
ing atoms is ~2n 2 J 2 . Hence, the contribution of spin i to 



the energy is Ei = — DlSf^Y" — 2n 2 J 2 S\ z ' • An incoming 

neutron typically causes a spin fli 
the energy of such a transition is 



(z) (z) 

neutron typically causes a spin flip 5- = 2 — > 5- = 1; 



AE 3D + 2n 2 J 2 , 



(2) 



FIG. 3: (a) Unit cell with seven occupied sites of which the 
central site is active and the six dimmed sites are held frozen 
in an antiferromagnetic configuration, (b) The same unit cell 
with all the occupied sites being active (for clarity, only two 
of them are marked with flipping arrows). The corresponding 
frozen environment is the next shell of neighbors (not shown 
in the figure). 

the edges of the unit cell (Fig. [3]). 

Following the same pattern, we model the disordered 
sample with the following Hamiltonian: 

(id) 

+ Ji ^ CiCjSi ■ Sj + J 3 Y e i e jSi-Sjll) 

((id)) 

where the Si are S = 2 spin operators that represent the 
Fe atoms, and where the substitutional disorder is rep- 
resented by statistically independent random variables ti 
which take the value one with probability x (the fraction 
of Fe atoms) and the value zero with probability 1—x (the 
fraction of nonmagnetic Zn atoms symbolized as empty 



and therefore to a first approximation the spectrum con- 
sists of nine evenly spaced zero-width peaks with a bi- 
nomial distribution of intensities. The anisotropy pa- 
rameter D determines the average position of the peaks, 
while the spacing between peaks is proportional to J 2 - 
Actually, a more accurate description can be obtained 
by averaging the immediate generalization of Eq. ^ 

AE w 3D - 2n x Ji + 2n 2 J 2 - 2n 3 J 3 , (3) 

over the respective number of neighboring sites < n\ < 
2, < n 2 < 8, < ri 3 < 4. In Fig. [4] we show the re- 
sult of this calculation with the parameters of the pure 
sample^ and (zero- width) intensities proportional to the 
products of the respective combinatorial weights. This 
figure shows how the coupling constants J\ and J 3 con- 
tribute to the effective spread of the peaks. In fact, com- 
parison with the fit to the experimental data suggests 
that all the pure sample parameters need to be modified 
to describe the highly diluted sample and, in particu- 
lar, that the anisotropy parameter D is too large. As 
we will see in the forthcoming discussion, the parame- 
ters that provide the best fit to the experimental data 
depend on the approximation scheme used to study the 
model Eq. (IT]). However, the semiclassical picture given 



4 




23456789 
E (meV) 

FIG. 4: Semiclassical approximation to the spectrum pro- 
vided by Eq. (|3| averaged over the respective number of neigh- 
boring sites with the pure sample parameters.^ The solid line 
is a fit to the experimental data of Fig. [I] corresponding to 
x = 0.31 and q = 0.5, and suggests that the pure sample pa- 
rameters need to be modified to describe the highly diluted 
sample. 



by Eq. ^ and illustrated in Fig. [1] remains qualitatively 
correct in the full simulation. 

Our approach to simulate the experimental spectra 
is a two-step procedure. In the first step we use the 
full Hamiltonian ([I]) but maintain the approximation 

Si ■ Sj s» s[ z ^ for the exchange terms, so that we can 
generate typical local environments by a classical Monte 
Carlo simulation. More concretely, we generate equilib- 
rium spin configurations at T = 10 K using a heat bath 
combined with a cluster method in lattices with 2 x 32 3 
sites with an Fe density x — 0.31. After we equilibrate 
ten such lattices (samples) we pick at random on each 
sample 1000 nonempty sites. We call this site, together 
with its n nearest neighbors, as illustrated in Fig.[3]jb), a 
dynamic shell, and the next shell of atoms [not shown in 
Fig. |3^b)], the local environment. Figure [3]( a) illustrates 
a simpler version of this idea, in which there is only one 
spin in the dynamic shell (the central spin marked by a 
flipping arrow) and there are six frozen atoms in an an- 
tiferromagnetic state that constitute the local environ- 
ment. 

In the second step of our procedure we use the full 
Hamiltonian ([T) with full quantum spin operators for 
each atom on the dynamic shell, while the nondynamic 
spins that constitute the local environment are kept fixed 
and act in effect as boundary conditions for the dynamic 
shell. Since the third component of the total spin for 
each dynamic shell 



(4) 



i dynamic 



commutes with the Hamiltonian, we find the ground state 
(G.S.) within the subspace corresponding to the Neel 
state S^ z ' — SnccIj which amounts to diagonalizing a 



square matrix with up to 3000 states, and the excited 
states by diagonalizing the Hamiltonian restricted to the 
subspaces SW = SnccI ± 1 allowed by the selection rules. 
The transition energies are simply the energy differences 



e. 



E, 



G.S. 



(5) 



The exact calculation of the intensities in this exper- 
imental setting involves the matrix elements of a rather 
complicated interaction Hamiltonian.^3 We settle for an 
estimate of the relative intensities of these transitions 
and consider the simplest possible interaction operator 
(in fact, one of the terms appearing in the full expres- 
sion), which is proportional to S± , where i denotes the 
central atom of the dynamic shell. The contribution of 
each transition ([5| to the total intensity is proportional 
to 



Kfc,£to =S NW ±l|£f ) |G.S.,S«= 1 S; 



Noel/ 



(6) 



We have found that the main effect of the simplified tran- 
sition matrix element ^ is to suppress the contributions 
of high-energy transitions. Note also that the combina- 
torial factor is already included in our sampling of the 
simulation results. Finally, we add up the contributions 
of all the possible transitions in our 1000 samples, cal- 
culate the convolution of this result with the measured 
instrumental resolution function (a Gaussian with mod- 
erate energy-dependent width), and normalize the result 
to match the maximum count of the experimental curves 
(it would seem better to match the integrated intensity in 
the experimental range, but as we will see, the resulting 
widths are too narrow). 

Figure [5] shows the results of these procedures for the 
two dynamic shells illustrated in Fig. [3] in the energy 
range between 2 meV and 8 meV. The figures show the 
experimental points, a numerical fitting to these points, 
and our simulation results. The one-site calculation re- 
produces quite well the main features of the experimen- 
tal spectrum, including the average position of the peaks 
(controlled by the anisotropy parameter D), although 
the separation between peaks (controlled by the cou- 
pling constant J 2) is too large and, as we anticipated, 
the widths are too narrow. 

Although the simulation results corresponding to the 
dynamic shell of Fig. [3^b) feature wider widths, the av- 
erage position is clearly shifted to high energies, which 
suggests that, not only the value of J2, but also the value 
of the anisotropy parameter D may be too large in this 
context. A possible explanation might be related to the 
method by which the pure sample parameters^ are ob- 
tained, whereby a semiclassical approximation is used to 
determine the spectrum parametrically as a function of 
D, Ji, J2 and J3, and later these parameters are fitted 
to match the experimental results. In essence, this pro- 
cedure involves a calculation to first order in 1/5, which 
should give better results for S > 1 i.e., for the one- 
site approximation. This possibility has already been 
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FIG. 5: Experimental data for q — 0.5, numerical fitting and 
simulated spectrum at T = 10 K with the parameters of the 
pure sample^ for: (a) a one-site calculation corresponding to 
the dynamic shell of Fig. |3|a); (b) calculation corresponding 
to the dynamic shell of Fig. [3jb) . 

noticed. For example, in Ref. [T3] certain empirical rela- 
tions between the parameters of the pure of the diluted 
sample are proposed. 

Following these ideas, in Fig. [6^ a) we show the result of 
an optimization of the parameters J 2 and D to match the 
experimental results, which yielded J 2 — 3.35 cm -1 and 
D = 5.25 cm -1 . Unfortunately, the correlation between 
these parameters and the uncertainties prevents a more 
accurate determination of these values or the simultane- 
ous optimization of the less significant J\ and J3 that we 
have kept fixed. Although the widths of the peaks are 
still too narrow, the intensities are quite well accounted 
for, and even a last peak at E sw 7.4 meV seems to be 
reproduced. 

Finally, as an estimate of the thermal effects in our 
simulations, in Fig. j6^b) we show a similar calculation at 
T = K, i.e., with a purely antiferromagnetic state (no 
thermal disorder) of the environment. Note the distinctly 
narrower widths and poorer intensity relations between 
the peaks, particularly at high energies. 

There are a variety of possible reasons to explain the 
larger experimental widths: our scattering operator is 
oversimplified, as it does not take into account the rel- 
ative orientation of the lattice and the wave vector of 
the incoming neutron, which we also assume perfectly 
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FIG. 6: Experimental points, numerical fitting and simulated 
spectrum with fitted values of D and J2 corresponding to the 
dynamic shell of Fig. ^b) for: (a) T = I K; (b) T = K 
(i.e., purely antiferromagnetic environment). 

well defined (i.e., we neglect the spread of the neutron 
beam); the dynamic shells and their environments have 
been obtained from a classical (rather than quantum) 
Monte Carlo, which surely overestimates the spin order- 
ing at low temperatures; and the dynamic shells are lim- 
ited to spins and their immediate neighbors. 



IV. DISCUSSION 

In summary, we have presented high-resolution spec- 
tra from neutron scattering experiments, conducted over 
Fe x Zn 1 _ x F 2 close to its percolation threshold. We model 
these spectra in terms of a site diluted Heisenberg model, 
containing both ferromagnetic and antiferromagnetic ex- 
change interactions. In spite of its simplicity and with 
only a moderate adjustment of the parameters, the pro- 
posed model accounts quite well for the position and in- 
tensity relations of the peaks in the spectra. This success 
is probably due to the validity of our main hypothesis, 
namely the local nature of the spin excitations in these 
systems which lie close to the percolation threshold for 
the Fe lattice. 

Whereas local spin excitations dominate the energy 
spectrum for x near x p in Fe!_ x Zn x F 2 , we cannot rule 
out a very small contribution from fracton excitations in 
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a similar energy range. Fracton excitations as well as lo- 
cal spin excitations coexist in isotropic systems and both 
may exist in the small anisotropy Mn 1 _ x Zn x F 2 system 
as x approaches x p . In that case, modeling the local 
spin excitations could aid in separating the two types of 
excitations, allowing the characterization of local spin ex- 
citations as well as the persistence of fracton excitations 
under conditions of weak anisotropy. 
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